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m-i, We calculate the corrections to the three-gluon decay widths T^ g of cc and bb 
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quarkonia due to the finite extension of the QQ — > 2>g vertex function. The 



K> . widths computed with zero range vertex are reduced by a factor 7 where 7 = 

SH : 

0.31 for the J /ip and 7 = 0.69 for the T. These large corrections necessitate 

a redetermination of the values a s (fj,) extracted from T^ g . We find a s (m c ) = 

0.28 ± .01 and a s (m b ) = 0.20 ± .01. 
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Since the discoveries of the charmonium and bottonium states, the physics of 
heavy quark systems has been of continuous interest |1], Q] . The bound state struc- 
ture of heavy quarkonia becomes simpler with increasing mass, because the con- 
stituents move more and more non-relativistically in a static potential which is 
dominated for r < 0.1 fm by the Coulomb-like field and for r > 0.2 fm by the confin- 
ing potential. The Schrodinger equation with this potential reproduces successfully 
the energies of the bound states ls,2p, .... The leptonic widths T(QQ — >■ /i/i) of 
the 3 S'i-states, modified by radiative corrections, can also be obtained by an expres- 
sion which is proportional to ^(O)) 2 , where ip(r) is the non-relativistic bound state 
wave function. The value |?/>(0)| 2 is rather sensitive to the potential || ||. Prob- 
lems exist with the gluonic width T(QQ — > 3g) and with the photo-gluonic width 
T(QQ — > ■ygg). One usually calculates these quantities in the approximation of a 
zero-range vertex, i.e. T(QQ — > 3g) oc |^(0)| 2 and analogous for T(QQ — * j2g). 
From a comparison with experiment one extracts a value for the strong interaction 
coupling constant a s (/x), where /i is the scale (fj, = m c for J/ip and jjl = mt, for T) 
[|, 0] . The analysis usually proceeds in three steps. In the first step, expressions for 
the ratios R\ and i?2 

T(QQ -» igg) 



Ri 



r(QQ ^ ggg) 



T(QQ -, 3g) 

are derived in leading order a s . In this order and using the zero- range approxima- 
tion for the vertex, the decay rates are obtained from analogous calculations for the 
three-photon decay of the triplet state of positronium with the photon-electron cou- 
pling replaced by the SU3 C coupling of gluons to quarks. Because of the zero-range 
approximation, the rate is proportional to |^(0)| 2 , the square of the non-relativistic 
wavefunction, which cancels out in the ratios R\ and i?2- In the second step, the lead- 



ing order result for R\ and R 2 is multiplied by the first order a s corrections, which 
represent a factor of the form (l + 5(/x)a s (/x)/7r). The coefficient B([x) depends both 
on the scale \x and the renormalization scheme. In a third step, one usually takes into 
account corrections of order (cf/m 2 ) phenomenologically ||. These corrections may 
come from the finite extension of the decay vertex and from relativistic corrections 
in the wavefunction or from other nonperturbative physics like condensates. In this 
paper we give the first ab initio calculation of the effect from the finite extension of 
the decay vertex, which we take into account by a calculated factor 7, where 7 = 1 
represents the zero range situation. The status of theoretical calculations for Ri 
and R 2 is given in Table 1 for the J/ijj and T together with the experimental values 
for the two ratios. Values for a s (n) extracted from a comparison of the theoretical 
expressions (7 = 1) with the data are also shown in this table. 

The following observations can be made for the zero-range case: (i) The ratios 
Ri and R2 give rather consistent values for a s (/x) for the same decaying system 
(J/ip or T). (ii) Using the values for R 2 which have small error bars, one obtains 
precise values for a s {m c ) and a s (rrib) from J/ip and T decays, respectively. The two 
values are nearly equal. This result is unexpected, because in a simple-minded QCD 
analysis, the scale dependence of a a (/i) would predict 

a s {m c ) ~ a s (m b ) J^|tT5 - °' 28 ± M ^ 

instead of a s (m c ) = 0.19 ± .01 where we used A = 200 MeV. This discrepancy has 
been observed before || and a correction factor has been parametrized but never 
calculated. 

We improve on the previous calculations by taking the quark-momentum de- 
pendence (or finite size) of the transition operator QQ — > 3g into account. This 
correction cancels in the ratio Ri, since the rates in the numerator and denominator 
are affected in the same way. In the ratio R 2 , only the three-gluon state is modi- 



fied, whereas the s-channel photon annihilation still samples |f/>(0)| 2 = | f d 3 qip(q)\ 2 , 
where ip(q) is the Fourier transform of ij){r). The three-gluon decay amplitude can 
be written as 

T(h, k 2 , h) = J j^jipBsiP, q)iM(p, q, k h k 2 , k 3 ), (3) 

where k iy i = 1,2,3 are the four momenta of the gluons, M is the expression for 
the QQ — > 3g vertex and ?Pbs(p,q) stands for the Bethe-Salpeter wavefunction of 
the quarkonium state. In eq. (§) p denotes the cm. momentum and q the relative 
momentum of the two quarks. In our calculation, we replace ipBS by the non- 
relativistic wavefunction ip(q) according to 

2n 
i>Bs(p, q) = —5(q°)u(p/2 + q)^{q)v{p/2 - q). (4) 

The reduction of the Bethe-Salpeter wavefunction to the nonrelativistic wavefunc- 
tion is a necessary approximation to proceed with the calculation. At short dis- 
tances, i.e. for high momenta \p\ > mq, it is expected that the relativistic wavefunc- 
tion falls off more slowly than the nonrelativistic wavefunction, since the relativistic 
transverse one-gluon exchange kernel is less damped. Only a part of this effect is 
taken into account in the radiative corrections || oc 0(a s ) which include transverse 
gluon exchange between the heavy Q and Q. The amplitude M factorizes into an 
expectation value for the colour part of the gluons with fixed colors i, j, k and the 
rest 

M(p,q,k u k 2 ,k 3 ) = Mijk Tr(^^^=\ 

= —j=(d ijk M s (p,q,k 1 ,k2,k 3 )+if ijk M a (p,q,k 1 ,k 2 ,k 3 ). (5) 

The first part M s includes the symmetric sum of the six diagrams arising from 
the interchange of the three gluons, while the antisymmetric part M a sums the 



same diagrams with a negative sign for the odd permutations, dijk and fijk are the 
structure constants of the group SU3. For q = 0, M a = 0, and one finds for the 
three-gluon rate the three-photon decay rate of positronium with the modification 
a 3 —* ^Q!g. For q 7^ 0, the contribution of M a is numerically found non-zero but 
negligible in comparison to M s . The vertex M is calculated from the expression 

M 123 (p, q, fa, fa, fa, ei, e 2 , e 3 ) = 

g^Ah ,,^L h h , /9 , , l , A-P fc (p,g)l. (6) 

1 - P/2+ /?+ Ai - m Q p/2+ jq- fa-m Q J 

Here, the vectors e^, i = 1, 3 describe the gluon polarization and P Sz (p, q) is the spin 
projection operator in the 3 Si state ||. We calculate the decay width T(QQ — > 3g) 
including the full ^-dependence in the matrix element M 

T(QQ - 3<?) = (7) 

1 , d3fc ld 3fcd3fc ^ ^ _ ^ 

pol 8u; i u; 2W3(27r) J 

-»-*-♦ /" d 3 q -»-»-» 

T(/ci, fc 2 , fca, ei, e*2, £3) = y -— ^M(j9, <f, fei, fc 2 , fca, ei, e* 2 , e* 3 ) ■ ^(g). (8) 

The extensive numerical calculation to determine T is performed in the following 
way. First we use the symbolic program REDUCE to evaluate the trace in eq. (§). 
The output of this symbolic program is a rather long function M (50 kbytes) of 
q, fa, fa, fa, ei, e* 2 , e* 3 . The T-matrix (§) is obtained after integrating M ■ ip(q) over 
Qx, %, Qz with an adaptive routine which gives a precision of 5%. We use wave- 
functions from a numerical solution of the Schrodinger equation with Richardson's 
potential || given in ref . Q . This potential goes over into the Coulomb potential for 
short distances and becomes a confining linear potential at large distances. Richard- 
son's potential gives a wavefunction ip(q) which falls off asymptotically as oi s {(f , )/\q\ 4: - 
The matrix element M(q) oc \q\ (cf. fig. 1) together with the wavefunction leads 
analytically to a very weakly diverging integral proportional to lnln|<f|rnax, which 



is in agreement with our numerical calculations. These show no sensitivity to a 
variation of the cutoff between |g|max = 3mQ and |g|max = 7mg. 

Finally the square of the T-matrix is summed over the gluon polarization vec- 



tors and gluon momenta fci = (0, 0, Ui), k 2 = (— yl — x 2 u 3 , 0, — uj\ — xlj 3 ), k 3 = 
(\/l — x 2 to 3 , 0,xu)3). Energy conservation in eq. (|7|) fixes x: 

X = - [(2m Q -Wi -uj 3 ) 2 -u\ -u)l). (9) 

The two polarization vectors of each gluon are chosen to be e*(l) = e y and e*j(2) = 
ki x e y /\ki\. We use a seven point two-dimensional integration M routine to do the 
final phase space integration over the Dalitz triangle given by < cu\ < m, < 
co>2 < m and u\ + 002 > m (cf. eq. ([7|)). We estimate our total numerical accuracy 
to be 10%. For consistency we checked our routines by setting q = in the matrix 
element M, which leads to the zero-range approximation T^°\ The corresponding 
decay width T°(QQ — > 3g) agrees with the well- known result for positronium decay 
modified by the gluonic couplings (eq. @) and the strong interaction wavefunction 
ip(r = 0). We denote by 

7 = T(QQ - 3g)/rl°\QQ - 3p) (10) 

the reduction factor for the three-gluon decay due to the finite size effect in the 
matrix element M(q). We obtain 

[ 0.31 ±0.03 for J/tp 
1 ~ \ 0.69 ±0.07 for T . ^ > 

At first sight this correction is surprisingly large, especially for charmonium. In 
order to demonstrate the origin of the finite-size reduction, we plot in fig. 1 the two 
factors in the integral eq. (§), the radial Gaussian wavefunction (fip(\q]) for J ftp 
and T as a function of \q\/m,Q where mq = m c for the J/ip and uiq = vn^ for the T, 
together with the matrix element M s (q, k%, /C2, k 3 , e*i, e*2, 6*3) which has been averaged 
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over the directions of q, but is still a function of the polarization directions. The 
matrix element M s (q) has its maximum at q = and falls off to a minimum at 
q/rriQ ~ 1. The exact shape depends on the gluon polarization vectors and the spin 
projection s z . We have checked the shape of the g-dependence of M s (q) for various 
combinations of momenta fcj and polarizations e*j. The shown case is rather typical. 
It is evident from fig. 1 that M s (q) cannot be well approximated by a parabolic 
function 

M s (q)=M s (0)(l-aq 2 /m 2 Q ) (12) 

over the whole domain of values q of interest. If we approximate ip{q) by a Gaussian, 
the form of eq. (|l^) would lead to a correction factor 

Y st « (l - a±^) , (13) 

where < v 2 > is the expectation value of the velocity of the quarkonium state under 
consideration |H| . The authors of ref . assume the correction factor to be of the 
form eq. flT3|) and fit the coefficient a. They find a between 2.9 and 3.5. For a = 3, 
the correction amounts to 7 = .28 for the J/ip and 7 = .78 for the T. These values 
are not too far from our numerical values eq. (|H]). This agreement is gratifying, 
since it suggests that the finite-size correction is probably the most important effect. 
The correction factor 7 calculated in this paper leads to new values of a s {m c ) and 
a s (mfc) if 7 is inserted into the theoretical expressions of R2 in Table 1. (The ratio 
Ri is unaffected by this correction.) The QCD corrections to the gluonic width for 
heavy mesons have been calculated in ref. || for the T-meson. Special care is taken 
in this reference to avoid double counting of the one gluon Coulomb exchange. The 
contribution to radiative corrections due to the Coulomb exchange between Q and 
Q must be dropped (i.e. class / in fig. 1 of ref. ||), because it is included already 
in the wavefunction. Only transverse gluons contribute to this class of radiative 



corrections. Ref. H provides a practical concise summary of all QCD corrections to 
the decays of quarkonia, also for cc-quarkonia. The values of a s (/i) deduced from the 
ratio R2 are shown in Table 1. The new values of a s calculated with 7 7^ 1 have to 
be compared with those from Ri and are consistent. Now the values for a s (m c ) and 
a s (mb) are rather different from each other, as is expected for the running coupling 
constant. We can use the values for a s (/i) to determine the QCD scale parameter 
A = A™-, which is related to a s {n) by 

a(u) = ^ (l 6(153 -IQn^lntln/zyA') \ 

slW (33-2n / )ln/i 2 /A 2 \ (33 - 2n/) 2 ln(/i 2 /A 2 ) ")' K } 

We use rif — 4 and \x = m c = 1.5 GeV for J ftp and fi = rrib = 4.9 GeV for T. We 

take the values of a s (fi) with the finite size corrections from Table 1 (for T we use 

an average of the values deduced from R\ and R2), insert them into eq. ( |I4"D and 

solve for A. The values of a s (/z) and the calculated scale parameters A are 

a s (m c ) = 0.28 ±.01 A§L = 200 ± 20 MeV , . 

K 5 (15) 

a s {m b ) = 0.195 ±.007 Ag^ = 255 ± 45 MeV. V ; 



We have also investigated the dependence on the choice of the scale \i discussed 
in ref. [llj by going from \x = tuq to \x = 2rriQ/3. The radiative corrections for 



the strong decay widths calculated from the expressions of ref. change rather 
drastically for the two cases. The new values of a s (2mQ/3) and the calculated 
values of the scale parameter Ai±. are the following 



,2 — ^ -0.35 ±0.01 A^=210±20 

T \ 

1„.\_noJinni A ( 4 ) 



u a (\m h ) =0.24 ±0.01 AgL = 300±40. 



Eqs. (^,|1^) summarize the main results of our paper. The values for A^± agree 



with the present world average [12| of A-^ = 238±30±60 MeV, where the first error 



MS 



is statistical and systematic and the second reflects the uncertainty in the scale. 

In this paper, we have shown that the inclusion of the finite size of the QQ —>■ 3g 
vertex reduces significantly the 3g decay widths of the J ftp and the T. The finite size 
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modifications necessitate a redetermination of the values for a s (fi). These new de- 
values are in agreement with the QCD evolution of the coupling constant a s (fi) with 
the scale \x and lead to consistent values of A. The corrections due to finite size of the 
3g vertex can be calculated without free parameter. Because of the large size of the 
reduction of the 3-gluon width of charmonium, the question arises whether there are 
not other relativistic corrections. Indeed we believe that the use of a nonrelativistic 
wavefunction is the most crucial approximation in our calculation. To improve on 
this approximation, much more conceptual and numerical effort would have to be 
spent both on the relativistic bound-state problem and the radiative corrections. 
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Figure Caption 

Fig. 1: The wavefunctions q 2 ip{q) for the Tq(i71q = nib, foil line) and for the 
J/ip{m,Q = m c , dashed line) are given in arbitrary units as functions of q/rriQ. 
The matrix elements M(q) (eq. ([5])) integrated over the angle of q are indicated 
for different spin directions of the 3 Si-quarkonia. Squares give M(q) for J z = — 1; 
crosses for J z — 1 and diamonds for J z = 0. The z-axis is defined by the momentum 
of one fast gluon. The energies and polarizations (u\ =002 = 0.97m,Q, u 3 = 0.06m<g 
and e*i = e*2 = £3 = (0, 1,0)) of the gluons are chosen to show typical variations of 
M(q) with q/rriQ. 
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Theory 


a s (n) 


j/4> 


0.10 ±0.04 


^2(1-3.02*) 


09 
19 ± 

.05 


T 


0.0274 ±0.0016 


#2(i_2.5^) 


0.189 ±0.011 
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p _ n»Si->ggg) 

«2 - rpsi-wO 


Exp 


Theory 


a s ()u);7 = 1 


7 


a s (n) 


w 


10.4 ±0.7 


5fe^(l + 1.59a./7r). 7 * 


0.190 ±.004 


0.31 ±0.03 


0.28 ±.01 


T 


32.5 ±0.9 


i°^K(i + .43a s /7r)- 7T 


0.181 ±.02 


0.69 ±0.07 


0.20 ±.01 



Table 1 : Experimental values and theoretical expressions for the ratios 
i?i and i?2 for the J ftp and T states. The experimental values are taken 
from Kwong et al. [5] for Ri(J/ip) and from the review talk by Kobel 
[7]. We used a~ l {m^) = 133.7 and a _1 (m M ) = 132.0 [6]. The coupling 
constants a s (/-0 are computed from the experimental ratios with the help 
of the theoretical expressions. Finite-range corrections cancel in the ratio 
Ri. In the case of i?2, the finite-range corrections are contained in the 
factors 7. The values for a s (fi) obtained in the zero-range approximation 
7 = 1 are compared with those which include the finite-range corrections. 



11 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9401233vl 



